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Abstract 

We consider the decoherence of phase space histories in quantum 
brownian motion models, consisting of a particle moving under a po- 
tential V{x) in contact with a heat bath of temperature T and dissipa- 
tion constant 7 in the Markovian regime. The evolution of the density 
operator for this system is thus described by a non-unitary master 
equation. The phase space histories are described by quasiprojectors 
consisting of gaussian density matrices smeared over large phase space 
cells. They are characterized by the size [r] of the phase space cell 
together with the size [M] of the margin (the region at the boundary 
of r in which the Weyl symbol of the projector goes from 1 to 0). By 
generalizing an earlier result of Hagedorn on the unitary evolution of 
coherent states, we show that an initial Gaussian density matrix re- 
mains approximately Gaussian under non-unitary time evolution, and 
derive a bound giving the validity of this approximation. This result 
is then used, following earlier work of Omnes to compute the time 
evolution of the phase space projectors under the master equation: 
The evolution of an initial projector P onto a cell T is approximately 
equal to another phase space projector P' onto a cell V , where V 
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is the classical dissipative evolution of T. Furthermore, the expected 
unpredictability due to thermal fluctuations is reflected in the fact 
that the margin of P' (and hence the effective region it occupies in 
phase space) is greater than that of P. We thus show that histo- 
ries of phase space samplings approximately decohere, and that the 
probabilities for these histories are peaked about classical dissipative 
evolution, but with an element of unpredictability due to the noise 
produced by the environment. 
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1 Introduction 



Ever since the first days of quantum mechanics, the question of the emergence 
of the classical deterministic world from the underlyinfg probabilistic laws is 
considered of particular importance. It is closely connected to the measure- 
ment problem, since any measurement involves the interaction of the system 
under observation with a large deterministic device. 

Different aspects of this problem have been adressed extensively in the 
literature. Still, there are issues that remain unclear. Among them are ques- 
tion of the validity of the description of many body systems with collective 
variables that evolve under quasiclassical equations of motion, the conditions 
under which classical predictability arises in a quantum system, its possible 
breakdown at long times and how it is affected by environmentally induced 
noise. In this paper we will address these problems, as appearing in a wide 
class of open quantum systems, namely quantum Brownian motion models. 

Our analysis lies within the framework of the decoherent histories ap- 
proach to quantum mechanics, as set out by Griffiths [0],Omnes 0, [| |j] and 
Gell-Mann and Hartle |J. The main element of this formalism is the notion 
of history. A history is a string of projection operators at consecutive mo- 
ments of time. We can build sets of histories, by taking a partition of the 
unit operators into projectors P a (t) at each moment of time. These operators 
satisfy: 



PctPfi — SafiPa 





a 



We usually consider them to evolve in the Heisenberg picture 



P a (t) = e lHt P a e 



-iHt 




An element of a set of such histories reads : 





A space of histories has a natural candidate probability measure 
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This measure can be interpreted as a probability only if the standard proba- 
bility sum rules are satisfied. The set of histories is then called consistent (or 
decoherent). The condition for this can be writen in terms of the decoherence 
functional , a complex valued function of pairs of histories: 

D(a,a>)=Tr(C aPo C a ,) (1.6) 

The necessary and sufficient condition for the sum rules to be satisfied is 
then 

ReD(a, a') = a ^ a' (1.7) 
We usually employ the stricter condition 

D(a,a') = a^a' (1.8) 

The terms consistency and decoherence of the set of histories are given in 
the literature to the conditions (1.7) and (1.8) respectively. That assignment 
of a probability measure to a set of histories, allows us to reason with the 
histories using the rules of classical logic. 

In general, decoherence appears in physical systems, only for coarse grained 
histories. A typical example of coarse graining is to ignore a number of de- 
grees of freedom of the system, or take projectors into large intervals of 
position or momentum. The decoherence is in almost all physical situations 
only approximate, that is, the probability sum rules are satisfied only within 
an order of e, where e is a small positive number. 

In quantum mechanics the most general proposition for a physical system 
at a given instant of time is represented by a projection operator. In classical 
mechanics the corresponding entity is the phase space cell. Thus, when 
we want to discuss classical deterministic behaviour, we have to construct 
entities that bridge the gap between the two formalisms. 

This is done through the introduction of quasiprojectors, that is, ap- 
proximate projectors onto phase space cells. They are constructed in such 
a way, that their corresponding symbol in the Wigner representation is a 
smoothened characteristic function onto some phase space cell. Omnes || 
has shown that histories corresponding to classical hamiltonian evolution 
of sufficiently large and regular phase space cells, approximately satisfy the 
consistency conditions. Also the conditional probability for our system to 
be found in the cell T 2 at time t 2 , if at time t\ it was within the cell r 1; is 
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almost unity, provided T2 is obtained from Ti through evolution according 
to the classical equations of motion. 

In this paper, we are interested in generalizing these results for a class of 
open systems, namely quantum Brownian motion models. These consist of a 
particle moving under a potential V(x) and coupled to a large system which 
is taken as the environment. Usually we take as an environment a collection 
of harmonic oscillators in a thermal state. By tracing out the environment 
degrees of freedom, we obtain an effective evolution equation for the density 
matrix of the particle. The coupling of the environment to the particle is 
contained within two effects: dissipation and diffusion. 

We concentrate on a particular case of those models, namely when the 
evolution equations are Markovian. This corresponds to the case of ohmic 
environment at high temperature. This is mainly due to the fact that only 
in the Markovian limit can we write the decoherence functional in terms of 
the reduced density matrix. 

These models can be used as prototypes for more complicated systems, as 
for instance particle detectors in early universe or continuous measurement. 
What is more interesting for our approach, is that we can model this way the 
effective behaviour of collective or hydrodynamic variables. We can identify 
the degrees of freedom of the particle with distinguished variables in a many- 
body system, like for instance the center of mass position and momentum of 
a macroscopic body, while the environment represents all the ignored degrees 
of freedom like the positions and momenta of its constituents. 

The issue we are concentrating on, is the conditions under which the 
description of the system with the classical dissipative equations of motion 
can be accurate. Essentially we want to prove approximate determinism for 
quantum dissipative systems. Our approach is rather different from other 
authors that have discussed decoherence of phase space histories || [7| or 
configuration space histories 

H and closer to the spirit of Omnes || || . 

Our main results can be summarized, as follows: We establish the degree 
of the validity of the gaussian approximation for the Markovian regime of 
quantum Brownian motion models. We show that a gaussian density matrix 
remains approximately gaussian centered around the classical path for a large 
class of potentials V(x). The validity of the approximation is quite good as 
long as the spread in position is much smaller than the scale on which the 
non-linearities in the potential become significant. 
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We use this result to show, that quasiprojectors onto sufficiently large 
and smooth phase space cells evolve according to the classical equations of 
motion. The environment induces a degradation on the evolved projectors, 
which leads essentially to a loss of predictability 

Finally, we construct the decoherence functional in the Markovian regime 
and establish that histories corresponding to classical evolution on phase 
space approximately satisfy the consistency condition, and thus approximate 
determinism arises. 

This paper is structured as follows: After giving a brief review of quan- 
tum Brownian motion models, we proceed to discuss phase space projectors 
in section 3. We generalize the construction of quasiprojectors from coherent 
states of Omnes, using gaussian density matrices. In section 4 we state a 
theorem on the degree of validity of the gaussian approximation in our mod- 
els, using which we establish how the quasiprojectors we constructed evolve 
in Section 5. We construct histories corresponding to classical evolution in 
Section 6 and discuss their decoherence properties. In the last section we 
summarize and discuss our results. 

The notation we will use is as follows: Whenever there is a possibility of 
confusion the reduced density matrix will be denoted as p, otherwise plainly 
as p. The Hilbert space of the wave functions will be denoted as H^ and 
the one of the density matrices under the Hilbert-Schmidt inner product as 
H p . We write the inner product in using the Dirac notation {ipi\ip2) and 
in Hp as (pi, pz) . The action of an operator C on an element of H p will be 
denoted as C[p\. The volume of a phase space cell T is written [r]. In Section 
4 and in the Appendix we will use units such that H is dimensionless. 



2 Quantum Brownian Motion Models 

We first give the basic features of the model, within which we are going to 
discuss phase space decoherence. We consider a particle of mass M moving 
under the influence of a potential V(x) and in contact with a heat bath of 
harmonic oscillators || [TO], [Ufl . 

The total action of the system reads : 

S tot [x{t),q n {t)} = f dt^Mx 2 - V(x)]+J2 I dt[\m n q n 2 
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-\ m nU 2 n q 2 n - C n q n x] (2.1) 

We concentrate in the behaviour of the distinguished particle and thus con- 
struct the reduced density matrix: 

am = n /«■**<*.*.«'.*) 

We can describe its time evolution using the reduced density matrix propa- 
gator J, defined by the relation: 

Pt(x, y) = f dx dy J (x, y, t\x , y , 0)p(x o , y ) 

Under the assumption that the initial density matrix po factorizes, the prop- 
agator is given by the path integral expression : 

J(x f ,y f ,t | x o ,y o ,0) = J DxDyexp(^S[x] - l -S[y] + ^W[x,y}) (2.2) 
where 

S[x] = f dt{-Mx 2 -V(x)) 
Jo 2 

W[x(t),y(t)] is the Feynman Vernon influence functional phase 

W[x(t),y(t)} = - fds [ S ds'[x( S )-y( S )}r]( S - S ')[x( S ')+y( S ')} 
Jo Jo 

+i f ds f ds'[x(s) - y(s)]u(s - s') [x(s') - y(s')} (2.3) 
Jo Jo 

The path integration is taken over all paths satisfying x(0) = xq, x(t) — Xf. 
The kernels rj(s) and v(s) are defined by : 

f°° du r . . , . Huj s , n JN 

u(s)= —I(u)coth(—)cosus 2.4 

JO 7T Ikl 



where 



V (8) = ^(s) (2.5) 
7 (a) = r^Monus (2.6) 

Jo 71 CO 



and I{uj) is the spectral density, 



ttC 2 

I{ U ) =Y,Ko> ~ <»n)~-f- (2.7) 



The kernel v(s) modifies the action of the distinguished system and leads 
to dissipation and renormalization of the potential. The kernel 77 (s) is respon- 
sible for noise and the process of decoherence. Both kernels are completely 
determined once a form for the spectral density is specified. A choice that 
proves to be convenient is 

I(u) = M^y- l exp(-^ 2 ) (2.8) 

where A is a cut-off, which is usually taken to be very large. 

The ohmic case (s = 1) is of particular interest. Taking A to infinity gives 

7(5) = M^5(s) 

which means that the corresponding classical equations of motion are local in 
time. The noise kernel is non-local for large A, except in the Fokker-Planck 
limit, kT » HA, in which we have, 

u(s) = —j—8(s) (2-9) 

In this regime the reduced density matrix p (we drop the tilde) satisfies 
the master equation 

dp = _ h ( d 2 p d 2 p dp dp 

dt 2MVdx 2 dy 2 ' 71 V)K dx dy> 

+ ^[V(x)-V(y)}p-^(x-y) 2 p (2.10) 
where we set D = 2M^ykT, or in operator form : 

= C[p] ' (2.11) 
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This master equation describes a Markovian process. The time evolution of 
the initial density matrix can be represented by the action of a one-parameter 
semigroup with generator C on the state space. 

We will represent an element of the semigroup as e ct and the time evolu- 
tion of an initial density matrix p as eCt [Po] ■ 
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3 Phase space projectors 



3.1 General considerations 

In order to discuss phase space histories, we have to find elements of the 
quantum mechanical formalism that correspond closely to the notion of a 
phase space cell, at least in the macroscopic domain. Since any proposition 
about a physical system at a given instant of time can be represented as a 
projection operator, what we need is a class of operators that project onto the 
subspace of physical states which have position and momentum well within 
some phase space cell T. 

Position and momentum cannot be simultaneously specified with arbi- 
trary accuracy, so we must consider approximate projection operators . We 
call an operator P a quasiprojector into a phase space cell T if any wave func- 
tion well localized inside T (outside T ) is an approximate eigenfunction of P 
with eigenvalue 1(0). The construction of such operators is particularly 
easy in the Weyl representation, in which we associate with P a function 
f(x,^) on a mock phase space, which we call its associated symbol 



In general, a function f(x,^) on a two dimensional phase space is called 
a symbol of order m, where m is real, if its derivatives of all orders have the 
following bounds 



where C a p are constants called the seminorms of / . Note also, that x and £ 
are in dimensionless units (scaled by some typical dimensions L and P ). 

One possibility is to take for the symbol corresponding to a quasi-projector 
on a phase space cell T the characteristic function of the cell T. This does not 
define a good operator because of the discontinuity at the boundary so in- 
stead, we consider a smoothened characteristic function [|J . That is, we take 
/ to be equal to 1 inside and to outside T, except in a small region along the 
boundary of T, where it interpolates smoothly from 1 to 0. This transition 
region is called the margin M of the cell T . In order for the operator to be 
well defined, the phase space cell has to be regular; loosely speaking it should 
not develope structure on the scale of h and the volume of the margin M 



f{ x ,0 = J dyexp(-i£y/h)(x + ||P|x - |) 



(3.1) 



<m^(2,0l < C a p(l +X 2 +^)(m-a-0)/2 



(3.2) 
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should be much smaller than the volume of V. The parameter e = [M]/[r] is 
a measure of the validity of the approximation. For a sufficiently regular cell 
it is of the order of (h/LP)? , where L, P are typical length and momentum 
scales for our cell. We will return to the 

issue of regularity in more detail later. The quasiprojectors P have a 
number of interesting properties ||: 

1. Their trace is proportional to the volume of the phase-space cell. 

TrP = [T]/(2nh) (3.3) 

2. They are very close to true projectors: 

Tr\P-P 2 \<ce (3.4) 

where c is a number of the order of unity. 

3. For two different projectors on the same phase space cell P and P' 
with corresponding parameters e and e' we have : 



Tr\P - P'\ < c(e + e') (3.5) 

4. If Ti, T 2 and Ti n T 2 are three phase space cells and P\, P 2 , P\2 
quasiprojectors associated to them, the operator SP = P1P2 — P12 is bounded 
in trace norm by : 

Tr\5P\ < K^\(h/LP) a (3.6) 
2nn 

where [r] = s«p( [IY], [r 2 ]) , a > ~ and K of order unity. 



3.2 Coherent state quasiprojectors 

One particular way to construct quasiprojectors is through the use of coherent 
states [Q. Consider the gaussian wave function : 

This induces a metric : 

*- p) = !* 2 + mW 2 (3 ' 8) 
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on phase space with respect to which one can give a precise definition of a 
regular phase space cell ||, |j. 

We call a phase space cell regular to order e, with respect to one partic- 
ular family of coherent states if: 

1. The curvature radii with respect to the metric d of dT are larger than 
I in absolute value. 

2. The margin of the cell M which is defined as : 

M= P| e{x,p,l) 

{x,p)edr 

where e(x,p,l) is the ellipsis defined by : 

e(x,p, I) = {(x' , p')\d(x — x',p — p') < I 2 } 

satisfies [M] < e[T) . 

3. The numbers e and 1 satisfy : e~ 21 < e. 

If T is a regular cell to order e then the operator : 
is a quasiprojector associated with T. We can readily see that its symbol : 

^« = / r ^ ex »t-|( 1+r2 )^-«) 2 

~;i ({ ~ p)2 ~l^ (I ~'' )(^ ~ ! ' ), (3 ' 10) 

is up to corrections of order e~ /2 equal to 1 (0) inside (respectively outside) 
T except for the margin M where it interpolates between those two values. 
It also satisfies the conditions for being a symbol of arbitrary negative order. 

3.3 Generalized gaussian projectors 

We are mainly interested in the time evolution of the quasiprojectors, when 
our system is under the influence of the thermal environment, as in the 
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case of quantum Brownian motion models. Due to the non-unitarity of the 
evolution, the operator in time t will take a form that cannot be related 
to the quasiprojectors (3.9) in a straightforward way, as in ||. Omnes has 
shown that under a large class of potentials, the density operators of the form 
\9gp) (9gp\ evolve approximately into operators of the same form but belonging 
to a different family of coherent states (different parameters S and r). The 
non - unitarity of the evolution in our case implies that an operator of this 
form become mixed. Therefore, we need a larger class of gaussian density 
operators. 

We propose a class of quasiprojectors, that are defined through gaussian 
density matrices, which we choose to parametrize by the set of five numbers 
(E,F,r,q,p) as 

l x \p\y) = (-?-) a exp[-4(^^ - q? 

~2~h} x ~ v ' ~^m^^i — v)( x -y) + i ]-( x -y)\ l 3 - 11 ) 

which is defined for S < AF in order to satisfy the positivity requirement. 
For these density matrices we can verify that, 

= (Trp 2 )! = (E/4F)* = (h/2A) < 1 (3.12) 
where A is the Wigner function area defined as : 

A 2 = (Aqf(Ap) 2 - Cl q 

with : 



C Pg = ^((q- (q))(p ~ <P» + (P - (p»(9 " (?»> 



1 

2 

A is a measure of the phase space area in which the gaussian density matrix 
is localized. 

This density matrix defines a metric on phase space in analogy to (3.8): 

d{X > P) = ^h x2+ 4F(ll§ri/ (3 ' 13) 

With respect to this metric, we can define the notion of the margin of a 
phase space cell and the notion of regularity (up to order e) as before . We 
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define the operator : 




dqdp 



p(E,F,r,q,p) 



(3.14) 



r 2irh 



This is a quasiprojector associated with the cell T. 
To see this, consider its symbol, 



/(*, = (E/4F)* / r ^ exp[-|(l + ^r 2 )(x - g) 2 




(3.15) 



It is a symbol of arbitrary negative order and has the form of the most 
general gaussian smeared characteristic function associated to the cell T. We 
can verify that it takes values 1(0) inside (outside) T up to corrections of 
order e~ l , where / is defined as before but with respect to the metric (3.13). 
It interpolates between those two values only in the margin, which has a 
width of order A^ . 

How can we compare this class of projectors to the ones defined through 
the general gaussian states (3.11)? A coherent state is localized within a 
volume of order h/2 and is unable to distinguish points found within a volume 
of this size. For this reason we expect the margin associated with a projector 
onto a phase space cell to have an area of order h and the parameter e to be 
of order (h/LP)z. This estimation is backed by a detailed calculation using 
microlocal analysis, that is, we can verify that this order of magnitude for the 
parameter e is the one that gives the best approximation of a quasiprojector 
to an exact projector |3|, [|. 

On the other hand the gaussian density matrix (3.11) is localized within a 
volume of order A. Therefore we expect the width of the margin to be of order 
A and the parameter e to be of order (A /LP) 2 for the class of projectors 
(3.14). This is clearly, not the "best" quasiprojector one can associate to a 
phase space cell, but this generalization is necessary. As long as the quantity 
(2A/h)z is of the order of unity we still have a good approximation for a 
phase space projector. The larger the value of A is, the worse the degree 
of approximation to a true projector is. We are therefore going to refer to 
projectors with minimum value of A = h/2 (the ones constructed from pure 
gaussians) as the maximum resolution projectors. 
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We will see that the time evolution of an initial generalized coherent state 
remains within a good approximation a gaussian density matrix, with time 
increasing value of A. We know from the study of linear systems that the 
value of A increases polynomially in the short time limit, to become constant 



for times much larger than the relaxation time 7 ||12|| . Its asymptotic value 
corresponds to the uncertainty due to thermal fluctuations, and is typically 
many orders of magnitude larger than the quantum ones. Still, for sufficiently 
large phase space cells A is much smaller than the size of the cell, the quantity 
e remains quite small and our operators (3.14) are still close to true projectors. 
On the other hand, for smaller cells for which h << LP < Aoo there comes 
some time (usually of the order of 7 _1 ), when the evolved quasiprojectors 
cannot anymore distinguish the phase space cell, since the size of the margin 
has become essentially as large as the whole of the cell. 



4 Time evolution of gaussian density matri- 
ces 

We are interested in the evolution of the gaussian density matrices (3.11) 
in the class of quantum Brownian motion models described by the master 
equation (2.11). As we mentioned, this expression is valid in the Markovian 
regime (high temperature and ohmic environment). 

For linear potentials, we know (see for instance || |12[ ) that the propa- 
gation is gaussian, and an initial gaussian density matrix remains gaussian 
centered around the classical path. Also in the case of zero coupling to the 
environment (unitary evolution) Hagedorn |13|] has established that a gen- 
eralized coherent state retains its shape for a period of time and its center 
follows the classical equations of motion. The error of this approximation is 
of the order of (h/LP) x , where A < |. 

We seek a generalization of these results. We want to show, that for a 
large class of potentials V(x) the master equation (2.11) respects the gaussian 
nature of the density matrices (3.11) and that their centers follow the classical 
equations of motion (with dissipation). We shall show that this is a good 
approximation as long as the spread in position is sufficiently small and that 
the error is of the order of (h/LP) x as in the purely hamiltonian case. 
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4.1 The theorem 

In this section, we present our result in the form of a theorem, a detailed 
proof of which can be found in the Appendix. We assume that the potential 
V(x) satisfies: 

1. V^ 2 \x) is continuous and uniformly Lipschitz on compact subsets of 
1Z (i.e. given any R> there exists k such that \V {2 \x) -V {2 \y)\ < f3\x-y\ 
whenever \x\ < R, \y\ < R. 

2. \V(x)\ <e Mx2 

3. V(x) is bounded from below. 

Assume that at t — the system is in the state po — p(^o, Fo, fo, Qo,Po)- 
The exact solution to the master equation (2.11) can be writen formally as : 

Pt = e°M 

where C is the generator of the one-parameter semigroup acting on the state 
space, which is determined by the dynamics. 
Then we have the following theorem : 

Theorem: For each T>0, iV>0, p < \, a <\—p and < t < T, 
such that 

(-^-r 1 < nk- 2p 

there exists C > 0, 5 > for which : 

\\e a [p ] -p(m,nt),r(t),q(t),p(t))\\HS < Ch x (4.1) 

whenever h < 5. We denote A = 3a — 1. 

The quantities (E(t),F(t),r(t),q(t),p(t)) are solutions to the system of 
differential equations : 

q = — (4.2) 

P = -V'(q)-2 1P (4.3) 

t = ±T?r (4.4) 

1 2D 
EFr - 4 7 F + — (4.5) 



M ' h 

Er 2 2 2 



r = 



F -2 1 r + -V {2 \q) (4.6) 



2M M E 
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under the initial conditions 



(E(0), F(0), r(0), g(0),p(0)) = (E , F , r , g , Po) (4.7) 

Note that we have used dimensionless units with L = 1 and P = 1. 

The proof is based on the observation that the generator L can be split 
into two parts, described in the Appendix: £ = S + U, where S is the 
generator that corresponds to the Brownian motion of a free particle and U 
a generator that contains only the effects of the potential. The action of S 
can be exactly determined (for instance using the density matrix propagator 
which can be exactly computed for the free particle), while the action of U 
can be approximated with a quadratic expression. Thus we get expressions 
for e [po] and e ut [po] which can be combined to give an expression for e ct [po] 
through the use of the Trotter product formula. The proof follows closely the 
one of Hagedorn |]13||, with the main differences being the non-unitarity of 



the generator S and that our Hilbert space is the one of the density matrix 
under the Hilbert-Schmidt norm. 

The proof involves taking a Taylor expansion of the potential and keeping 
only the quadratic terms. Therefore our results are exact for the harmonic 
oscillator and the free particle case. 

4.2 Uncertainties 

We can write the uncertainties in position Aq and in momentum Ap as well 
as the correlation C pq associated with the gaussian density matrix in terms 
of the parameters E, r and F: 

(Aq) 2 = h/E (4.8) 
[Apf = hF(l + ^r 2 ) (4.9) 

Tit 

C m = -- (4.10) 

From this we can see the interpretation of the parameters E, F and r we used 
to parametrize the gaussian density matrices. Clearly E _1 is a measure of the 
spread in position and r of the correlation between position and momentum. 
E/F is proportional to Trp 2 and therefore is a measure of the non-purity 
of the state. The original parameters E, F, r have been more convenient in 
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the course of the proof, but in order to have a clearer interpretation of the 
results, we shift our attention to the set of variables Aq, Ap, C pq . Within 
our approximation they evolve according to the equations: 



d_ 

di 



(A,) 2 = fa 



pq 



(4.11) 



It 



(Ap) 2 = -4 7 (Ap) 2 - 2C pq V i2 \q) + 2D 



(4.12) 



d 
dt 



M 



1 



(Ap) 2 -2 1 C pq -(Aq) 2 V^(q) 



(4.13) 



We notice, that the diffusion coefficient appears only in the equation 
for the momentum uncertainty, and at short times is the dominant term. 
This means that the spread in momentum is more effective than the corre- 
sponding spread in position, which depends on the diffusion coefficient only 
indirectly. It is also important to stress the reluctance of particles with large 
mass to undergo a rapid growth of the fluctuations. 

4.3 Validity of the approximation 

The essential requirement for the validity of our approximation is that the 
quantity remains sufficiently large through the evolution according to 

the equations (4.2-4.6). Since we have that £ < AF, this is equivalent to 
saying that S remains sufficiently large (much larger than Nh 2p ). But a 
large value of £ means that the particle is well localized in position. We can 
thus say, that while the state remains well localized in position (that means 
Aq < N-2h~ p ), the gaussian approximation will be quite good. When, 
mainly due to the diffusion, the state of the particle has become spread 
in space, the large scale structure of the potential becomes important and 
even weak non-linearities will contribute significantly in the evolution, thus 
rendering the gaussian approximation invalid. Note, that the constants A 
and p satisfy : 



which shows that the larger our tolerance for the spread in position, the 
larger the error stemming from our approximation. 
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In general, we cannot say much about the time when the gaussian ap- 
proximation breaks down, since this depends crucially on the form of the 
potential, the particle mass and the position and spread of the initial gaus- 
sian. 

Excluding cases of potentials varying significantly in the microscopic scale 
or potentials where tunnelling effects can be important, we expect the gaus- 
sian approximation to hold within a very good accuracy for times much 
smaller than the typical time scale in our dynamics : 7 _1 . For larger times, 
the validity of the approximation, depends heavily on the length scale L on 
which V( 2 '(x) varies. If this length scale is much larger than the size of the 
thermal fluctuations the approximation will hold for timescales some orders 
of magnitude larger than 7 _1 . If this is not the case, the approximation will 
have broken down much earlier. Within a few times 7 _1 the particle will 
tend towards the state of thermal equilibrium. 

We can clarify those ideas by examining the simplest case of a system 
moving in a potential exhibiting only weak nonlinearities : 



This effect of the non-linearities becomes significant at length scales of 
the order of 



Assuming weak coupling to the environment, the uncertainty in position is 
given in leading order to e and 7/0; : 



with (Aq') 2 containing the effects of the unitary evolution. The gaussian 
approximation will break down when Aq ~ L. We can readily verify that for 
7] « u A /kT, Aq will remain much smaller than L, while for larger values of 
77 the approximation will break down at a time scale of the order of 7 -1 . 

We should note, that for a large variety of physical systems, it is not 
necessary to assume weak nonlinearities, in order to have a large value of 
L. Systems with potentials corresponding to a spatial average of many mi- 
croscopic degrees of freedom have typical values for L that can be said to 
correspond to a macroscopic scale. 




(4.14) 




77 




(4.15) 
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A relevant question is to give an estimation of the error term Ch x of our 
approximation. For the case of the closed system, a very rough estimation is 

h x f dt\V i3) (q(t))\Aq(t) 



o 

in dimensionless units 0. When considering the open system case we can see 
from the detailed study of the proof, that a term of Trp 2 < 1 enters the right 
hand side of the inequalities, thus rendering our approximation better. Also 
for U-shaped potential the effect of the dissipation is to make the values of 
the third derivatives of the potential smaller in absolute value. For the time, 
that the increase in the uncertainty due to the diffusion is small, we expect 
our approximation to be better when taking into account the coupling to the 
environment. On the other hand, at times larger than the typical time where 
thermal fluctuations overcome the quantum ones, the more effective spread 
of the density matrix makes the effect of non-linearities more important, 
and thus reduces the accuracy of the gaussian approximation. Therefore, we 
arrive at the following picture: 

At times less than (h/jhT)^ the gaussian approximation is better when 
the environment is taken into account. 

At larger times because of the effects of diffusion, the approximation 
becomes gradually worse, since diffusion begins to affect strongly the uncer- 
tainty in the position. 

Finally we should make a remark on a problem that may arise when 
having a closer look at the proof of the theorem. We know that the master 



equation (2.11) does not preserve positivity at very short time scales [12, 14 
This is a result of taking the cut-off A to infinity. Our proof, being based 
on a discretization of the time and taking the continuous limit in the end, 
might be inadequate for a class of initial states. We can avoid this problem, 
by adding a term r][p, [p, p]} in the master equation (like the one appearing 
in [15]). The master equation thus becomes positive and this term does not 
change the nature of the proof. In the end we can set rj equal to zero, an 
approximation that is valid in the high temperature regime. 
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5 Evolution of quasiprojectors 



In this section we study the time evolution of the quasiprojectors P intro- 
duced in Section 3. We shall show that as long as the corresponding phase 
space cells remain large and regular, they evolve according to the classical 
equations of motion with dissipation. One expects that the coupling to an 
environment will induce a degradation in the quasiprojectors as they evolve. 
This is reflected in an increase of the size of the margin, which implies a loss 
of predictability. 

We work in the equivalent of the Heisenberg picture for open systems. 
That is, we assume that the density matrix is time independent and that the 
effect of the evolution is contained in the operators. Let us further assume 
that in this picture a (bounded) operator P evolves under the action of an 
one parameter semigroup with generator M.. The simple correspondence 
between Ai and C that exists in unitary evolution is lost here. 

We determine Ai by demanding that the probabilities are the same in 
both pictures. This means : 

Tr(e Mt [P]p) = Tr(Pe Ct [p}) (5.1) 

which is equivalent to demanding : 

Tr(M[P]p) =Tr(PC[p}) (5.2) 

Having the expression (2.11) for C it is easy to show that : 

dP 1 v 2 

- jjr{P,[P,x\}-§l*,lx,P}\ (5-3) 

The relation of this equation to the one for p is more clarifying in the 
position representation : 

~W(x) - V(y)]P(x, y) - 2 7 P(x, y) - ^(x - y) 2 P(x, y) (5.4) 
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Comparing this equation to the one for the p, we notice the change of 
sign in the terms corresponding to the dynamical evolution and the appear- 
ance of the extra term —2^P. Because of this term the equation fails to 
preserve the trace. This is an expected feature, when studying the evolution 
of quasiprojectors, since classical dissipative evolution does not preserve the 
phase space area. 

We can study now the evolution of gaussian operators W(E, F, r, q,p) : 

It is straightforward to show, using the treatment found in the appendix, 
that within an error of C(h/LP) x in the Hilbert-Schmidt norm, 

e Mt [W(Z ,F ,r ,q ,p )} = J(t) W(E(t),F(t),r(t),q(t),p(t)) (5.6) 

where E(t), r(t), F(t), q(t), p(t) are the solutions to the equations, 

(5.7) 
(5.8) 
(5.9) 

(5.10) 

(5.11) 

under the conditions 

(£(0), F(0), r(0), g(0), p(0)) = (E , F , r , q , p ) (5.12) 

The equations for q and p are the backwards classical equations of 

motion (not the time reverse) and J(t) is the Jacobian of the transforma- 
tion from (q ,p ) to (q(t),p(t)) 

d{q(t),p{t)) =e2lt 
d(q ,Po) 



Q = 


p 
M 




p = 


2 7 p + V\q) 




£ = 


L r 

M 




F = 


— — EFr + 47F + 
M 


2D 

~T 


r = 


£r 2 2 „ 

1 F + 2<yr • 

2M M 1 


~^\q) 
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Now let us consider the evolution of a quasiprojector Pp, associated with 
a phase space cell T which is regular to order e. Consider an operator that 
at t — has maximum resolution in phase space. We thus have at t — 

Pr = /r^T ^( S o,£ /4,r ,go,Po) (5-14) 
As it evolves under the action of e Mt this projector becomes 

e Mt [P T ] = Jit) W(E(t),F(t),r(t), qd (-t),p d (-t)) (5.15) 
within an error of order 

(h/LP) x TrP r 

Performing a transformation from (qo,po) to {qd{—t),Pci(—i)) we readily ver- 
ify that : 

e Mt [Pr] = ^ ^wMt),F(t),r(t),q,p) = P Tt + 0(e') (5.16) 

where T t is the phase space cell obtained from T through the classical equa- 
tions of motion and e' = [A/ LP) 2 determines the degree of regularity of T t . 
Here we have: 

is the Wigner function area. P' is clearly another quasiprojector 
associated with the phase space cell T t . Essentially we have : 

He^fPr] - P'Mhs < c(h/LP) l iTrP r (5.18) 

From (3.4) we know that within an error of e', Pf is equal to Pp t , where Pp t 
is the quasiprojector defined through the gaussian operator with parameters 
(S ,4S ,r ) 

||^r t - PtMhs < Tr\P Vt - P Tt \ < de'TrP Tt (5.19) 

Assuming that during the evolution the error of the gaussian approxima- 
tion is less than e', we conclude from (5.18) and (5.19) that 

||e^[P r ] - P Tt \\ HS < Ce'TrP T (5.20) 
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This is our main result: A projector Pp onto some phase space cell T 
evolves under the action of the dynamical semigroup into a projector Pp t as- 
sociated with the phase space cell T t , which is obtained from T through the 
classical dissipative equations of motion. This is an approximate result, ac- 
curate within an error of order e'. This parameter increases with time solely 
because of the noise induced by the environment. This implies that the mar- 
gin of the quasiprojector increases with time evolution and the predictability 
becomes gradually worse. 

We can say that the effective region in phase space occupied by a pro- 
jector Pr consists of the cell T , plus the margin. Thus the effective region 
the projector occupies increases in a way that has nothing to do with the 
corresponding classical equations of motion ( unlike the shrinking due to dis- 
sipation which is a classical feature). A measure for this increase, which can 
be attributed solely to the diffusion, is 

P [r t ] + [M] l + e K LP J V W ' ' V ; 

We can view this result as follows: Sufficiently coarse grained observables 
(projectors onto large phase space cells) evolve under classical deterministic 
equations plus the action of noise terms. These terms, stemming from the 
ignored degrees of freedom of both the system and the environment induce an 
error of order e' to the classical equations of motions. The parameter measur- 
ing the size of the noise terms and correspondingly the loss of predictability 
is the Wigner function area A ||12|| , which in general is an increasing function 
of time. 

It may be shown that in the short time limit (t << 7 _1 ) for all initial 
operators (5.14) it is given by (T! 



^Hy + f 1 !^ 4 (5-22) 

We obtain the same result, when considering the backward time evolution 
of the quasiprojectors. 



6 Histories 
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6.1 Construction of the decoherence functional 



Using our previous results we now construct the phase space histories for the 
Brownian particle and study their decoherence properties. 

We have found how the quasiprojectors associated with phase space cells 
evolve in time. We can use this knowledge in the calculation of the decoher- 
ence functional. Our results are of use only when it can be written as: 

D(a, a') = Tr s \P n a K\ : _ x [. . . J* Kg [p]P^ . . .]P£J (6.1) 
in terms of a superoperator K\ in the space of the reduced density matrices. 



This is actually valid in the Markovian regime [15]. 

In the Caldeira -Leggett model the environment is taken to be essentially 
infinite and thus the environmental degrees of freedom remain very close to 
the state of thermal equilibrium. The correlations created will in general be 
of the order of A -1 and since the evolution has no memory their effect will 
be negligible in a time scale of order A -1 . Remembering that in our regime 
we have assumed kT >> HA, we conclude that if the time interval between 
the two projections satisfies : 

t 2 -t 1 > 0(H/kT) 

we can write the decoherence functional for two-time histories as 

D(a, a') = S^TriP^e^-^P^fcPbJ (6.2) 

or using the definition of the generator M. 

D(a,a>) = Sa^Trie-W^P^P^falfy (6.3) 

We can assign a probability measure 

p(C a ) = Tr((e- M ^\Pl\P^p tl P^) (6.4) 

to the history C a if the probability sum rules are satisfied. A sufficient con- 
dition for this is the vanishing of the off-diagonal elements of the decoherence 
functional. 
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6.2 Approximate decoherence and determinism 

The sum rules are not satisfied exactly in our case, but only within an error 
of order e. A good condition for this approximate consistency is : 

|rr( e -*(*-*>[Pi]Pi/*Pi,)| < 

6 £Tr( e -*<*-*>[J*]Pi/^i) ( 6 - 5 ) 

Let us concentrate on the particular case of histories that correspond 
to the evolution of phase space cells along the classical trajectories. That 
means, we consider the set of histories C\ = (Pi,P 2 ) and C 2 = (Pi, P2) 
where Pi and P2 are quasiprojectors associated with the cells I\ and T 2 and 
T 2 is the cell obtained from the classical evolution of IV In order to achieve 
minimum error in our estimations we restrict our quasiprojectors to have 
minimum resolution in phase space. The condition (6.8) reads : 

e (Tr((e- M ^- tl )[P 2 ]P 1 p tl Pi)+Tr(e- M ( t - tl )[P 2 ]P 1 p tl P 1 )) (6.6) 

From the analysis in the previous section, we know that within an error 
of order (A(t 2 - h) / 'LP)\TrP x 

e- M{t ^[P 2 }= P 1 (6.7) 

This means that the left hand side in (6.6) is of order 

(A(t 2 -t 1 )/LP)^TrP 1 

and the right hand side of order 

eTr(Pip tl ) < eTrP 1 

Thus the probability sum rules are satisfied within an order 

e = (A(t 2 -t 1 )/LP)* 

A takes values between h/2 and its asymptotic value which is of the order of 
the thermal fluctuations in phase space. Therefore, for a macroscopic phase 
space cell e is a small number and decoherence is good. It is clear that the 
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degree of decoherence is better in the case of a hamiltonian system. We can 
thus compute the conditional probability, that the particle being within I\ 
at time t 1 will be within T 2 at time t 2 . This is clearly, 

p(T 1 ,t 1 ^T 2 ,t 2 ) = l-0(e) (6.8) 

Therefore up to an order of e there is good agreement between the predictions 
of classical and quantum dynamics for the evolution of the Brownian particle. 
This agreement is best when the phase space cells are regular and the time 
evolution preserves this regularity. These results are independent of the 
initial state of the system. One might expect that the degree of decoherence 
would be in inverse relation to the degree of predictability. Here, they are 
essentially the same, as in the case where the environment is not present. 

The extension of this result to n-time histories is straightforward. Con- 
sider the history C = (Pi, t±] P 2 , t 2 ; . . . ; P n , t n ) where P; is the maximum 

resolution quasiprojector associated with a cell Tj and r i+ i is obtained 
from Tj 

through the classical equations of motion. We construct from C a set of 
histories by replacing one or 

more of the Pi 's with Pj, and define the parameter e = max(A(ti — 
ti-i)/LP)2. If we assume that for the intermeddiate times 

min(t m - U) > 0(k/kT) 

and that the dynamics preserve the regularity of the cells to an order e, then 
this set of histories satisfies the probability summation rules to this order. 

7 Conclusions 

In this paper, we have studied histories corresponding to the evolution o f 
phase space cells in quantum Brownian motion models with ohmic dissipation 
and Markovian dynamics. We constructed these histories using the general- 
ization (3.14) of the quasiprojectors used by Omnes. The estimation (4.1) 
for the validity of the gaussian approximation has enabled us to establish, 
how those quasiprojectors evolve. 

We showed that as long as the phase space cell is large and regular, his- 
tories corresponding to the evolution according to the classical deterministic 
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equations of motion approximately satisfy the sum rules. The order of mag- 
nitude of the corresponding error depends only on the Wigner function area 
A (a function only of the time difference between projections). A(t) can be 
computed as the solution of the system (5.7 - 5.11) under the initial condition 
(5.12). It starts as containing purely quantum uncertainties 

at t — and at a timescale of 7 _1 becomes of the order of magnitude of 
its asymptotic value, which corresponds to thermal fluctuations. 

There are three timescales in our model, which determine qualitative 
changes in the description of the dynamics: the Markovian time tM = h/kT, 
the decoherence time = (h/^kT)^ and the relaxation time t r = 7" 1 . For 
times t << tu dynamics are not Markovian and therefore we cannot talk 
about predictability. For times t such that < t < t d , sufficiently large 
phase space cells evolve almost under deterministic equations of motion and 
the details of the potential are not of importance in the increase 

of the fluctuations. At t — td the thermal fluctuations overcome the 
quantum ones \T2\. For t > the degree of decoherence and predictability 



gradually worsens and is destroyed for cells with area 

less than the thermal uncertainty. For these cells any sense of predictabil- 
ity has been lost at t — t r . In the contrary sufficiently macroscopic cells, 
continue to evolve within a good approximation according to the classical 
equations of motion even at times larger than r y~ 1 . Eventually at long times, 
approximately deterministic evolution will break down, because the dissipa- 
tive nature of the evolution tends to shrink the phase space cells, until their 
area is of the order of magnitude of the thermal fluctuations. We therefore see 
that the breakdown of predictability is dependent mainly on the parameters 
of the environment and the structure of the potential is largely irrelevant, 
unlike the unitary case. This comes probably from the fact that we have 
considered an one dimensional problem. It is to be expected that in a sys- 
tem with more degrees of freedom, time parameters of the potential (i.e the 
Lyapunoff exponents of the classical solutions) will play a more important 
role in the determination of when classical predictability breaks down. 

Perhaps contrary to our expectations, the coupling to the environment 
tends to make the decoherence properties of histories worse. This means that 
histories corresponding to evolution of phase space cells are not the ones that 
give the sharpest correspondence with a deterministic classical description. 
This is to be expected, since the classical description of quantum brownian 
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motion is that of a stochastic process in phase space. 

Essentially, there are two mechanisms that produce decoherence. One is 
the interaction with an environment and the other is the existence of phase 
space projection operators that have an intrinsic almost deterministic time 
evolution. In the case of purely hamiltonian dynamics only the latter ap- 
pears. In our case, we have both, but it is clear that again the latter is 
the dominant one. This is attributed to the simplicity of the correspond- 
ing classical equations, in which the effect of the environment is contained 
only within the term — 2 r yp. It is therefore natural that observables close 
to the ones used in purely Hamiltonian systems will decohere giving rise to 
approximate predictability. 

With those results we have achieved three things: First, we have general- 
ized the results of Omnes on classical determinism for a class of open systems. 
In fact, we can generalize our results are valid for any Markovian equation of 



the Lindblad type [18], with environment operators L n linear in position and 
momentum. Second, we have verified the validity of classical equations of 
motion with linear dissipation, as a consequance of the underlying quantum 
mechanical evolution and sufficient coarse-graining. Third, taking quantum 
Brownian motion as a toy model, we have obtained a picture of the behaviour 
of collective variables in many body systems and the quasiclassical behaviour 
of these variables. 

Finally, we would like to compare our results with the ones obtained in 
17]] using the quantum state diffusion picture for the quantum Brownian 



motion. For the case of two-time phase space histories C = (ri,^) and 
under the assumption that they decohere, the authors gave for the associated 
probability the expression : 

p(C) = dp 2 dq 2 dpidqxJ(p2,q 2 ,t 2 \pi,qi,t 1 )f(p 1 ,q 1 ,t 1 ) (7.1) 

where f(p,q,t) is a classical probability distribution satisfying the Fokker- 
Planck equation and J is the associated propagator. We can verify that this 
result is equivalent to ours in the particular case of an harmonic potential, 
where the propagator is gaussian centered around the classical path. If the 
volume of the phase space cells is much larger than the Wigner function area 
A(t) associated 



with the propagator |12| , then the propagation is within an error of order 



(A/ LP) 2 a 5-function around the classical path and the conditional 
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probability is found 

p(ri,*i->r 2 ,f 2 ) = i -0(e) 
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A Proof of the theorem 

We are going to give the proof for the theorem stated in section 4. It follows 
the reasoning of the similar proof of Hagedorn for the case 

of closed system, where the reader can refer for further details. 

The evolution equation for the density matrices in the Markowian regime 
(2.11) corresponds to the action of an one-parameter semigroup with gener- 
ator C on the state space . We can write : 



where S and U are generators of one parameter semigroups given by : 



Splitting £ into S and U will enable us to use the Trotter product formula 
to determine e ct in terms of e st and e ut . 

We want to establish, how the class of states described by the gaussian 
density matrices (3.11) evolve under S and U. 
1. Evolution under S 



c = s + u 



%] = ^-[p 2 /2M,p] + l[x, {p,p}} - ^[x, [x,p\\ 
U [p] = hv{x),p] 



(A.l) 



(A.2) 
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Let us consider the solution of : 
with initial conditions 

p = p(S ,Fo,r ,go,Po) 

at time t — 0. This is just the solution of the master equation for the free 
particle. For short times t (t < r « 7~ x ) we have that : 

|| e 5 *[p(So,F ,r ,go,Po)] - p(Z,F,r,q,p)\\ H . s . < K{t/rf (A.3) 

where : 

q = qo + ^t (A.4) 

P = Po - ZlPot (A.5) 

£ = Eo + (^Egr )t (A.6) 

1 2D 
F = F o + (^S F r - 4 7 F + — )t (A.7) 

r = r ° " ~ M F ° " 27ro)t (A ' 8) 

2. Evolution under U - quadratic potential 

Consider the case of the most general quadratic potential : 

V{x) = a{x-q)+ l -{x-qf (A.9) 

We can easily verify : 

e"*[p(Eo,Fo,r , q ,Po)] = p(So, F , r + bt, q ,p - at) (A.10) 



3. Evolution under U - general potential 

For the case of a general potential satisfying the conditions we stated in 
section 4, we define : 

W q (x) = V(q) + V\q) + l -V {2 \q){x - qf (A.ll) 
32 



Denote by Uy and Uw q the generators corresponding to V(x) and W q (x) 
respectively. We have : 

\\e^[p] - e^[ P ]\\ 2 H.s < %{P, {U v -U Wq f[p\) (A.12) 

in the Hilbert-Schmidt inner product. Writing : O = \Uy — Uw q \ we see that 
the right hand side reads : 

t 2 2t 2 

¥ (0[p],0[p}) = ^ Tr(p 2 2 - pOpO) (A.13) 

where O is the operator \V(x) — W q {x)\ acting on the Hilbert space H^. In 
the coordinate representation this quantity reads : 

2t 2 f 

J dxdyp(x,y)p(y,x)0(x)[0(x) - 0(y)] < 
2t 2 f 

— s- / dx(x\p 2 \x)0 2 (x) = 

h J 

2t 2 S . nh . i r , 1 EF x 2n ^2/ \ /* ,n 

-^r — H^rr; / dx exp [ — — =-7- -(x - q) }0 2 (x) (A.14) 

The hypothesis that V(x) s uniformly Lifschitz implies that for qeK (some 
compact subset of 1Z) and \x — y\ < 1 there exists (3 such that: 

\V(x)-W q (x)\ <P\x-y\ 3 

Following the treatment of Hagedorn we can split the integral into two sec- 
tors: D = {x; \x — q\ < Ti a } and D = 1Z — D where a < \. 
The first part reads: 

2t 2 £ , irh , 1 r . 1 £F 



^27r/T£/4 + F' Vj? ' ft£/4 + F 

2t 2 „„ E , 7r/i s 1 /• . r 1 £F 

2t 2 



< ^(Trp 2 ,)^ 6 " (A.15) 
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The second part gives : 

2t 2 E , nh .1 r . r 1 SF 



y s ^ exph ^E74TF (a; - g)2]o2(a;) - 



h 2 27ra v S/4 + F y J D — -^ ftE/4 + F 1 
2t 2 E , Trft s i /■ , r 1 EF 



ft 2 1 2irh K E/A + F J Jd Fl ftS/4 + F v 



2t 2 E , Trft ,i /■ , rr 1 EF n71/rl/ , 21 

-5-C? — -( — )3 / drrexpj — -— h2M x-o 2 } 

ft 2 1 27rn v E/4 + F > ' is P1[ 2ft E/4 + F JV H! 1 



1 EF 

B ^ eXph 2lE74T^ ( ^ g)2] (A ' 16) 



where in the last step we used the Cauchy-Schwarz inequality at L 2 (D) (the 
space of real- valued square integrable functions on D). This quantity is 
smaller than 

(2Trp 2 )5 exp[Mft 2a - ft 2 "" 1 ] (A. 17) 



ft 

Assume that 



E F 

1 <Nh~ 2p (A. 18) 



V E + 4F 

for some > and p < a. This implies the existence of 5 > such that: 

E F 1 

M - ( j/r 1 < M — h^/N 2 

V E + 4F ; ~ 7 

for ft < <5 and taking a < | — p there exists constant C2 such that : 

exp[Mft 2a - h^+^/N 2 } < c 2 k 6a 

Thus we establish that when ( s+fp ) -1 — Nh~ 2p , qeK, for each AT > , 
p < \, a, < \ — p and t > there exists C > and 5 > , such that for 
ft < 5 : 

\\e Uvt \p\ - e Uw ^\p] \\ 2 H _s < Ch 3a ~\Trp 2 )h (A. 19) 

We now have expressions for the evolution induced by the generators S 
and U. We will combine them by using the Trotter product formula. 
4. Conclusion of the proof 

The hypotheses on V(x) ensure the existence of a bounded solution to 
the equations (4.2-4.6) for t <T. We can rewrite our previous result as : 

\\ e Uvt [po] _ eUw^p^ K C , h X /3T (A>2Q) 
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where A = 3a — 1 < \. We discretize the system (4.2-4.7) taking time step 
St = t/N . We have : 

QN(n) = q + ± Pj j^-5t (A.21) 

n 

p N (n) = po-E[^WW) + 2 7PJvW]^ (A.22) 
i=i 

Z N (n) = Z + ±f;[Y, 2 N (i)r N (i)}5t (A.23) 

i=l 

n 1 2D 

F N (n) = F + ^[—i: N (i)F N (i)r n (i)-A 1 F N (i) + —]8t (A.24) 

" r S w («H( 2 ) 2F N {i) o ... 

2 \/( 2 )(g^«)]5t (A.25) 



(A.26) 

Due to uniform convergence, we can always find Ni such that for N > Ni. 

\\p(X N (N), F N (N),r N (N), q N (N),p N (N)) 
-p(X(t),F(t),r(t),q(t), P (t))\\ HS < C'h x /3 (A.27) 

The Trotter product formula ensures the existence of N 2 such that for 
N > N 2 : 

\\(e Ct - [e St l N e ut ' N \ N )[p \\\ HS < C'h x /3 (A.28) 
For N > max(Ni, N 2 ) we have : 

\\ [e st/N e ut/N _ e st/N e w qo t/N M ^ Fo , ro , qo ,p )]\\ HS < C'h x /3 (A.29) 

We also have : 

\\ e st/N e ut/N [p ^ Fo> ro> ?QjPo)] _ p(I!jv ( 1 ) > ^(1), r^l), q N (l),p N (l))\\HS 

< \\[e St / N e ut / N -e^e^l^p^F^q^Us 
+ ||e' s '/ Ar p(So,F ,r + V^(q )t/N, q ,p - V'(q )t/N) 

-p(E N (l),F N (l),r N (l),q N (l),p N (l))\\ HS 

< C'h x /(3N) + Kit/rfN- 2 < Ch x /(3N) (A.30) 
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for N > N 3 > max(N u N 2 ) and with C > C 

For these values of N, we can obtain by iteration : 

n 

II II e st/N e m/N [p(J: , F , r , q , Po )} p(S iV (n), F N (n),r N (n), q N (n),p N (n) \\ HS 



< nCh x /(3N) (A.31) 



Taking the above inequality for n = N and using (9.27) and (9.30) we arrive 
at: 

||e £ * [p(S , F , r , q ,Po)} - HQ, r(t), q(t),p(t))\\ HS < Ch x (A.32) 

QED 
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